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Abstract — We develop a general downlink model for multi- 
antenna heterogeneous cellular networks (HetNets), where base 
stations (BSs) across tiers may differ in terms of transmit 
power, target signal-to-interference-ratio (SIR), deployment den- 
sity, number of transmit antennas and the type of multi-antenna 
transmission. In particular, we consider and compare space 
division multiple access (SDMA), single user beamforming (SU- 
BF), and baseline single antenna transmission. For this general 
model, the main contributions are (i) ordering results for both 
coverage probability and per user rate in closed-form for any 
BS distribution for the three considered techniques, using novel 
tools from stochastic orders, (ii) upper bounds on the coverage 
probability assuming a Poisson BS distribution, and (iii) a 
comparison of the area spectral efficiency (ASE). The analysis 
concretely demonstrates, for example, that for a given total 
number of transmit antennas in the network, it is preferable 
to spread them across many single-antenna BSs vs. fewer multi- 
antenna BSs. Another observation is that SU-BF provides higher 
coverage and per user data rate than SDMA, but SDMA is in 
some cases better in terms of ASE. 



I. Introduction 

Cellular networks are undergoing a significant makeover 
to handle increasing data demands prompted by the rising 
popularity of data hungry applications, such as video stream- 
ing and real time video calls J2). Handling these changing 
usage trends requires much higher data rates than contem- 
porary cellular networks were designed for. Two strategies 
that stand out to provide such spectral efficiencies are i) 
deploy low power nodes to reduce frequency reuse distance, 
and ii) equip BSs with multiple antennas to enable the use 
of multiple antenna techniques, such as beamforming and 
SDMA. Multiple antenna techniques are already relatively 
mature, being part of multiple wireless standards such as IEEE 
802.1 le WiMAX and 3GPP LTE-A (3j, apart from plethora 
of theoretical research activities in academia (4). Similarly, 
the concept of deploying low power nodes, also termed as 
small cells or HetNets, has been researched both in industry 
and academia for a fairly long time, see for example (5), 
[6 1 and the references therein. The standardization activities 
for HetNets have also started in 3GPP release 10 (7). These 
activities clearly indicate that multi-antenna techniques and 
HetNets will coexist and complement each other in the future 
wireless networks and should not be studied in isolation, as has 
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been typically done in the literature. In this paper, we address 
this problem and develop a general tractable model and the 
corresponding analytical tools for multi-antenna HetNets using 
techniques from stochastic orders and stochastic geometry. 

A. Related Work 

A potentially undesirable effect of small cells deployments 
is the increasing randomness in the BS locations due to lack of 
site planning. The departure of cellular networks from regular 
deployments towards nearly random deployments means that 
the cellular models that were based on the regularity assump- 
tion of BS locations, e.g., deterministic grid model, are no 
longer applicable. A more natural approach to model HetNets 
is by using random spatial models, where the locations of the 
BSs are assumed to form a realization of some point process, 
the simplest being the Poisson Point Process (PPP) (8), |9). 
This model has the advantages of being scalable to multi- 
ple classes of overlaid BSs and accurate to model location 
randomness, especially that of the small cells. Additionally, 
powerful tools from stochastic geometry can be used to derive 
performance results for general multi-tier networks in closed 
form, which was not even possible for single-tier networks 
using conventional models ||9). While sufficient progress has 
been made in modeling single-antenna HetNets (9j-(T2j, the 
efforts to understand multi-antenna HetNets have just begun, 
e.g., see fl3) . 

The main challenge in modeling multi-antenna HetNets is 
the number of possible multi-antenna techniques to choose 
from in each tier along with their tractable characterization. 
As a result, most prior works on multi-antenna HetNets have 
focused only on two-tier networks. For this paper, the most 
relevant one is fl4| , where SU-BF was shown to achieve 
better coverage than multiuser linear beamforming on the 
downlink of femtocell-aided cellular network assuming perfect 
channel state information (CSI). Random orthogonal beam- 
forming with max-rate scheduling and coordinated beamform- 
ing for femtocell underlay networks was analyzed in [15|, 
p6| , respectively. The effect of channel uncertainty on linear 
beamforming in two-tier networks was investigated in fT7) . In 
addition to the contributions in cellular networks, there has 
been extensive work on analyzing multi-antenna techniques 
in wireless ad hoc networks, which is also related to our 
work since several tools and techniques developed therein 
can be employed and extended to HetNets. Several single- 
user MIMO techniques, such as spatial diversity, open loop 
transmission and spatial multiplexing, have been studied, see 
for instance p8)-p0). The performance of multiuser MIMO 
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communication in a Poisson field of interferers, with perfect 
and quantized CSI at the transmitter was investigated in pT[ 
and J22) , respectively. 

B. Contributions and Outcomes 

Downlink model for multi-antenna HetNets: In Section II, 
we develop a comprehensive downlink model, which captures 
the heterogeneity of current deployments in terms of transmit 
powers, deployment densities and target SIRs across BSs of 
different tiers. It additionally models specific characteristics of 
multi-antenna HetNets, such as different number of transmit 
antennas across tiers, different transmission techniques and 
hence different number of users served across tiers in each 
resource block. We also consider the possibility of closed 
subscriber group or closed access in which a typical user is 
granted access to only a few BSs, while the rest purely act as 
interferers. 

Ordering results for coverage and rate: For general system 
models, such as the one considered in this paper, it is not 
always possible to express key performance metrics such as 
coverage probability and per user rate in closed form. In 
the absence of simple analytical expressions, it is difficult to 
compare different transmission techniques in general HetNet 
settings. To facilitate this comparison, in Section III we derive 
ordering results for both the coverage probability and the rate 
per user in both open and closed access networks, using novel 
tools from stochastic orders. Interested readers can refer to 
p3) , [24 1 for application of stochastic orders to conventional 
wireless networks. While circumventing the need for deriving 
coverage and rate expressions, it leads to several system design 
guidelines, such as i) to serve the same number of users as 
the number of antennas in each resource block, deploying 
single-antenna BSs provides higher coverage and per user 
rate than SDMA due to higher proximity gain, and ii) with 
more antennas than the users in each resource block, SU- 
BF provides even higher coverage and per user rate due to 
the additional beamforming gain. The BS locations for this 
analysis may be drawn from any general stationary point 
process, not necessarily independent across tiers, which is a 
significant generalization of earlier HetNet models based on 
the PPP assumption. 

Coverage probability upper bound and ASE comparison: 
While comparison of different configurations of multi-antenna 
HetNets in terms of coverage probability and average rate per 
user is conclusive from the ordering results, it does not directly 
capture the fact that some transmission techniques serve more 
users than the others and hence provide higher sum data rate. 
In order to capture this effect, in Section IV we additionally 
consider ASE, which gives the number of bits transmitted 
per unit area per unit time per unit bandwidth. To facilitate 
the comparison of transmission techniques in terms of ASE, 
we first derive an upper bound on the coverage probability 
of a typical user in both open and closed access networks 
and show that it can be reduced to a closed form expression 
for the "full" SDMA case (where the number of users served 
is equal to the number of antennas). For this derivation, we 
assume that the BS locations for each tier are drawn from 



independent PPPs. The tightness of the bound is studied and 
it is shown that the closed form bound derived for full SDMA 
is tight down to very low target SIRs. Using this expression 
and the previously known coverage probability expression for 
single-antenna HetNets, we derive ASE results for various 
transmission techniques in closed form. Main consequences 
of this analysis are: i) for the same density of BSs, single- 
antenna HetNets have lower ASE than SDMA since they 
serve fewer users. Interestingly, despite serving fewer users, 
SU-BF outperforms SDMA in moderate and high target SIR 
regime, and ii) when the BS densities are adjusted such that all 
the transmission techniques serve the same density of users, 
the ASEs of SU-BF, baseline single-antenna transmission 
and SDMA follow the same ordering as that of coverage 
probability and average rate per user. 

II. System Model 

A. System Setup and BS Location Model 

We consider K different classes or tiers of BSs, indexed 
by the set JC = {1,2,..., K}. The BSs across tiers differ in 
terms of their transmit power Pk with which they transmit 
to each user, deployment density A^, target SIR [3k, number 
of antennas Mk and number of users served by each BS in 
a given resource block \I>j. < M/., Vfc 6 JC. While in open 
access networks a typical mobile can connect to any BS, in 
closed access networks the access is restricted to B C JC tiers. 
For the ordering part, the locations of BSs of each tier are 
drawn from a general stationary point process The point 
processes $j. are not necessarily independent. Those familiar 
with the recent advances in modeling cellular networks with 
random spatial models will immediately recognize that this is a 
significant generalization over the known models based on the 
PPP assumption. This is enabled by the fact that the ordering 
results are based on the ordering of channel distributions 
for various setups and do not depend upon the spatial point 
process governing the locations of the BSs. However, we do 
require further assumptions for the ASE comparison, which 
involves the derivation of explicit expressions for coverage 
probabilities. Therefore, for the ASE analysis we will consider 
the more familiar independent PPP model, where each tier 
of BSs is modeled by an independent homogeneous PPP of 
density A^. 

Nevertheless, this is a fairly general model that captures 
the current deployment trends in 4G networks, e.g., it is 
easy to imagine hundreds of femtocells coexisting in each 
macrocell, transmitting at orders of magnitude lower power 
than macrocells, having relatively small number of antennas 
due to smaller form factor, serving smaller number of users 
due to smaller coverage footprints and providing restricted 
access to their own users due to a smaller backhaul capacity 
or privacy concerns. A two-tier illustration of the proposed 
system model is shown in Fig. [T] where a high power macro 
tier with four transmit antennas per BS coexists with a low 
power pico tier with two antennas per BS. Owing to its bigger 
coverage footprint, each macro BS serves higher load than 
its pico counterpart (25). The coverage region of each BS in 
this illustration corresponds to the region where it provides 
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Fig. 1. An illustration of a possible two-tier multi antenna HetNet configu- 
ration, with four antenna macro BSs serving two users and two antenna pico 
BSs serving one user each. 



the maximum average received power, thereby leading to a 
weighted Voronoi tessellation j9). 

We model the user locations by an independent PPP $ u and 
focus on the downlink analysis performed at a typical single- 
antenna user located at the origin. For the interference, we con- 
sider full-buffer model where the interfering BSs are assumed 
to be always transmitting (SJ, (9j. More sophisticated load 
models J25) along with non-uniform user distributions p6[ 
can also be considered but are out of the scope of this paper. 



• SDMA: In this case, a k tier BS with M k antennas 
serves ty k > 1 users in each resource block. When ty/. = 
M k , we term it as full SDMA, for which A fc = 1. 

• SU-BF: In this case, a k th tier BS serves ^ = 1 users 
in each resource block. 

• Single Input Single Output (SISO): Baseline single- 
antenna case |9), where each BS serves one user in each 
resource block. 

For each transmission technique, the received power at the 
typical single-antenna user located at origin from the BS 
located at x k E &k is 

Pr = Pkhk Xh \\x k \\- a , (1) 

where a is the path loss exponent, because we assume a per 
user power constraint in this formulation. The received SIR 
can now be expressed as: 



J2keK.^2y£<S>j\x k Pj9jy\\y\\ a 

For notational simplicity, we assume that the thermal noise is 
negligible as compared to the self interference and is hence ig- 
nored. This is justified in the current wireless networks, which 
are typically interference limited |28) . As will be evident from 
our analysis, thermal noise can be included in the proposed 
framework with little extra work. For cell association, we 
assume that the set of the candidate serving BSs is the 
collection of the BSs that provide the strongest instantaneous 
received power from each tier to which a typical mobile is 
allowed to connect. A typical user is said to be in coverage if 
the received SIR from one of these candidate serving BSs is 
more than the respective target SIR, as discussed in detail in 
the next section. 



B. Channel Model 

Before going into the technical details, it is important to 
understand that the channel power distribution of the link from 
a multi-antenna BS to a typical single-antenna user depends 
upon the transmission technique and whether it is a serving 
BS or an interferes For example, if it is a serving BS, it 
may precode its signals for a typical user depending upon the 
transmission technique, which may lead to a different channel 
distribution from the case when it simply acts as an interferes 
Therefore, to develop a general framework in which the BSs 
across tiers may differ in terms of the number of transmit 
antennas and the transmission technique, we assume that the 
channel power for the direct link from a k th tier BS located at 
Xk € R 2 to a typical user located at origin is denoted by h kXk 
and for the interfering link from a j th tier BS located at y g R 2 
is denoted by gj y . In this paper, we focus on zero-forcing 
precoding, under which for Rayleigh fading it can be argued 
that the channel power distributions of both the direct and 
the interfering links follow the Gamma distribution ]27| . As 
discussed in detail in Appendix |Aj it can therefore be shown 
that for zero forcing h kx ~ r(Aj,, 1) and gj y ~ T(^j,l), 
where Afc = M k — ^ + 1. For concreteness, we will focus 
on the following three transmission techniques in this paper: 



III. Ordering Results for Coverage and Rate 

This is the first main technical section of this paper, where 
we compare the performance of various transmission tech- 
niques in terms of coverage probability and average rate per 
user. We first study coverage probability in detail and then 
show that the analysis can be easily extended to study rate per 
user. We begin by formally defining the coverage probability. 

Definition 1 (Coverage probability). A typical user in an open 
access network is said to be in coverage if its downlink SIR 
from at least one of the BSs is higher than the corresponding 
target. This can be mathematically expressed as: 

P c = P ( I J max SIROfc) >(3 k ) . (3) 

The coverage probability can be equivalently defined as the 
average area in coverage or the average fraction of users in 
coverage. For closed access networks, the definition remains 
the same, except that the union is now over the set of tiers 
B C K, to which a typical user is allowed to connect. 

Remark 1 (Open access vs. closed access coverage). The 
coverage probability in open access networks is always higher 
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than in closed access networks. It follows directly by definition 
of coverage probability along with the fact that for B C /C 



ill max SIR(x fe ) > /3 fc 



vfces 



< 1 



I J max SIR(x fe ) > /? fe 



(4) 



\keic 



Owing to the complexity of the system model considered 
in this paper, it is not always possible to express coverage 
probability in simple closed form for any general system 
configuration. As evident from our analysis in the next section, 
this presents the first main challenge in comparing various 
transmission techniques. In this section, we take a slightly 
different view of this problem and focus on "ordering" the 
performance of different system configurations using novel 
tools from stochastic ordering. This powerful approach allows 
insights into the relative performance of different transmis- 
sion techniques, while circumventing the need to evaluate 
complicated expressions for the performance metrics such as 



coverage and rate. Interested readers can refer to |29| for 
details on stochastic orders. We begin by defining first order 
stochastic dominance as follows. 

Definition 2 (First order stochastic dominance). For any two 

random variables (rvs) Z\ and Z-i, Z\ (first order) stochasti- 
cally dominates Z% if and only if 



P[Zi > z] > ¥[Z 2 > z], Vz. 



(5) 



Equivalently, Z\ is greater than Zi in the usual stochastic 
order and is denoted by Z\ > st Z% 

Therefore, Z\ > st Z 2 if and only if the complementary 
cumulative distribution function (CCDF) of Z\ dominates 
that of Z 2 over the whole range. It is intuitively clear and 
will be made precise later in this section that the proper 
understanding of the ordering of received SIR for different 
system configurations plays a central role in studying their 
coverage and rate ordering. Now note that the received SIR 
can be alternatively expressed as: 



SIR(xfe) = 



Pk\\Xk\ 



p j9jy/h kx k \\y\\ 



(6) 



where the randomness due to propagation channel is lumped 
into ratios of Gamma random variables hk Xk and gj y . There- 
fore, it turns out that it is important first to understand the 
ordering of the ratios of Gamma random variables, which is 
studied next. 

A. Ordering of the Ratios of Gamma Random Variables 

For concreteness, define the ratio of two random variables 
Xi ~ T(k, 1) and X 2 ~ T(m, 1) by Z k , m = X 1 /X 2 . It is 
easy to derive the distribution of Z k m using basic algebraic 
manipulations and is given by the following Lemma. The proof 
is skipped. 

Lemma 1 (Distribution of the ratio of Gamma rvs). The 
cumulative distribution function of Zk. m — X\ / X2, where 



X\ ~ r(fc, 1) and X 2 ~ T(m, 1) is 

1 r(m + i) 



F z 



1 



1 



E 



r(m) ^ r(i-M) tz + i) m + i ' 

i—0 



(7) 



The form of this distribution function is such that for a given 
k\, mi and ki, fn%, it is not easy to derive conditions on these 
variables under which the CCDF of one ratio Zk,m dominates 
that of the other over the whole range of z. Therefore, the 
above result is of little help in providing direct ordering of 
two random variables Z kl _ mi and Zk 2 ,m 2 . We take an indirect 
route, which uses the following technical result about the 
equivalence in distribution of the Gamma random variable and 
the sum of i.i.d. exponentially distributed random variables. 
Given this technical result, the main result can then be proved 
using coupling arguments. We first state the equivalence result, 
which can be easily verified using characteristic functions. We 
then remark on how to use coupling arguments to establish 
stochastic dominance before stating the main result. 

Lemma 2. For i.i.d {Xi}, with Xi ~ exp(l), the random 
variable X = ^7J i=1 Xi is X ~ T(fc, 1). 

Remark 2 (Using coupling to establish stochastic dominance). 

One way to prove Z\ > st Z2 is to find two random variables 
Z\ and Z2 with the same distributions as Z\ and Z2, 
respectively, such that it is always the case that Z\ > Z^. This 
approach of using the same source of randomness to generate 
two random variables Z\ and Z\ satisfying the above relation 
and thereby establishing the stochastic dominance result is 
termed as coupling [30]. 



We now prove the following result on the ordering of the 
ratios of the Gamma random variables. 

Lemma 3 (Ordering of the ratios of Gamma rvs). A random 
variable Zk 1 .m 1 defined as the ratio of two Gamma random 
variables (first order) stochastically dominates Zk 2y m 2 if ki > 
h,2 and m\ < m,2- 

Proof: Using the equivalence in distribution of the 
Gamma random variable and the sum of exponential random 
variables shown in Lemma [2] we can generate a random 
variable Zl m with the same distribution as Zk i mi as 
follows 

i=l y M 



J k\ ,mi 



(8) 



j=l Y 2,i 

where {Y m n } is the set of i.i.d. random variables such that 
Y m ,n ~ exp(l). This equivalent representation will facilitate 
the use of standard coupling arguments, under which the goal 
now is to generate another random variable Zt _ with the 

e> K 2 .m 2 

same sources of randomness as that of Z^ , which has the 
same distribution as Zk 2 , m2 , and show that Z% m > Z^ 2 . 
Under the condition ki > fc 2 , this can be achieved by 
expressing m as follows 



J k\ ,ra\ 



> 



Ek 2 v x^ kl V 

j=l *2,i 

sr^k 2 v 
3=1 Y2 ' { 



(9) 



(10) 
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access multi-antenna HetNets. As evident from the analysis 
and remarked later in this section, all the results and insights 
carry over to closed access networks as well. 

Theorem 1 (Coverage ordering in open access networks). The 
coverage probability of a K-tier open access HetNet with 
system parameters {A&} and {^k} is higher than or equal to 
the one with system parameters {A' fe } and {^' k } if A& > A' fe 
and < W k for k € K. 

Proof: By definition, the coverage probability for open 
access networks can be expressed as 

P <= E1 | U max ^H^O^ > ft] 



(14) 



Fig. 2. The CCDFs of the ratios of Gamma random variables for different 
shape parameters k and m. 
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where (a) follows from the condition mi < m2, and (6) from 
Lemma [2] This completes the proof. ■ 

Remark 3. The set of conditions k\ > k^ and m\ < mi is 
stronger than the single condition k\jm\ > kilvni, which 
first comes to mind from the equivalence of Gamma random 
variables and the sum of exponential random variables stated 
in Lemma^ In fact it is easy to argue that the above stochastic 
dominance is not always true if the only condition on the 
variables is fci/mi > kilm%. For instance, consider a case 
where &2 S> k\ and m,2 S> mi such that k\jm\ > kijmi. 
The distribution of Zk 2 ,m 2 ' s concentrated around its mean 
and cannot be dominated by the distribution of Zk 1>mi due 
to significant difference in their shapes. This is illustrated in 
Fig. [2] where the ratio of Gamma random variables with shape 
parameters k\ = 4 and mi = 2 does not dominate the one 
with k<i = 100 and m,2 = 100 due to concentration, although 
kx/mi > kilrni holds. 

To prove the main ordering results of this section, we need 
to extend the stochastic dominance result of two random 
variables to multi-variate function of random variables. The 
result is given in the following Lemma and follows directly 
from the coupling arguments [30 1 . 

Lemma 4. If Xi > s t Yi for all i, then 



El (J 
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k£K y£$j\x k 



(15) 



where with slight abuse of notation (dropping the BS location 
from the subscript), Zjk = gjy/hk Xk is defined as the ratio 
of two Gamma random variables corresponding to the K-tier 
HetNet with system parameters {Aj.} and {v&fe}. Denote the 
corresponding ratio for the other system setup by Z'j k . By 
Lemma [3] Z' jk > st Z jk if ^ > * fc and A^ < A fc . Now 
note that the indicator function inside the expectation in ( [13] ) 
is element wise decreasing function of Zj k . Therefore, by 
Lemma [4] the proof follows. ■ 
Using this result, we can make some general comments 
about the coverage probability in certain realistic deployments. 
We begin by studying the effect of the number of users served 
in each tier on coverage probability. 

Corollary 1 (Effect of number of users). For two different K- 
tier open access HetNets, with the same number of antennas in 
each corresponding tier, the one that serves less users in each 
tier than the other provides higher coverage due to higher 
beamforming gain. 

The proof of the above corollary directly follows from the 
fact that under the same number of antennas for two setups, 
the one that serves less users in each tier has A& > A' fc 
and "Jfc < ^>' k for each tier, leading to higher coverage. An 
important extension of the above corollary is the comparison of 
the SDMA with SU-BF systems when the number of transmit 
antennas in each corresponding tier are the same. The result 
is stated as the following corollary. 



E[ ff ({JQ})] > E[g({Fj)] 



(13) 



for all multi-variate functions g 
each component. 



Corollary 2 (SDMA vs. SU-BF). For two K-tier HetNets, 
with the same number of antennas in each corresponding tier, 
that are non-decreasing in one performing SU-BF in each tier and the other performing 

SDMA, the coverage probability in the SU-BF case will always 
be higher. 



B. Coverage Probability Ordering 

Using Lemma [4] we now derive the following general 
result on the coverage probability ordering in K-tier open 



Another consequence of this general ordering result is the 
comparison of SISO with SU-BF and SDMA, with the SDMA 
case specialized to full SDMA. Recall that in case of full 
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SDMA, A fc = 1 and * fe = M k for all the tiers. The result is 
given in the following corollary. 

Corollary 3 (SU-BF vs. SISO vs. full SDMA). For three K- 
tier HetNet setups, one performing SU-BF in each tier, another 
doing single-antenna transmission in each tier and the last 
one doing full SDMA in each tier, the coverage probability 
in case of SU-BF is higher than that of SISO, which in 
turn outperforms full SDMA. The number of antennas in the 
corresponding tiers of SU-BF and full SDMA HetNets need 
not be the same. 

The proof follows from the fact that the shape parameters 
in case of SU-BF are A k = M k and = 1, where Mk > 1 
is the number of antennas; in case of SISO are A' k = ^>' k = 1; 
and in case of full SDMA are A% = 1 and ^' k = M k , where 
Mk > 1 is the number of transmit antennas. 

For closed access networks, it can be shown that the 
coverage ordering result derived in Theorem [T] holds under 
a slightly weaker condition because a typical mobile is not 
allowed to connect to all the tiers. The result is given as the 
following Corollary of Theorem [T] and the proof if skipped. 
Due to the similarity of this result with the open access case, 
the insights gained for the open access networks above carry 
over to the closed access networks as well. 

Corollary 4 (Coverage ordering in closed access). For two 
HetNets, with B C K, open access tiers, the one with system 
parameters {A^} and {^k} has & higher or equal coverage 
than the one with system parameters {A' fc } and {^j.} if A k > 
A' k for ke B and W fe < % for k € K.. 

C. Ordering Result for Rate per User 

Another metric of interest for the performance evaluation of 
HetNets is the rate achievable per user. In addition to the link 
quality (characterized in terms of SIR), it also depends upon 
the effective resources allocated to each user. For tractability, 
we make following two assumptions on resource allocation: i) 
each k th tier BS serves same number of users, and ii) each BS 
allocates equal time-frequency resources to all its users. For 
SDMA, it should be noted that several users will be scheduled 
on the same time-frequency resource block. Interested readers 
can refer to |3"T| for more details on the motivation and 
validation of these assumptions. Under these assumptions, we 
denote the effective time-frequency resources, e.g., bandwidth, 
allocated to a user served by a k th tier BS by Ok- The two 
assumptions on resource allocation ensure that Ok is the same 
for all the users served by any k th tier BS. Therefore, the 
downlink rate of a typical user served by a k th tier BS located 
at Xk € can be expressed as: 



R(x k ) = O k log 2 (l + SIR(xk)). 



(16) 



Due to the difficulties in modeling exact load on each BS [9|, 
which often requires characterization of the service areas for 
different types of BSs, it is challenging to characterize Ok 
and hence derive exact expressions for per user rate distribu- 
tion 1 3 1 1 . However, we now show that to compare different 
multi-antenna transmission techniques, this characterization is 
not required and the general ordering result derived above for 



the coverage probability can be easily extended for the rate 
per user as well. Before going into the technical details, it 
is important to note that the loading across tiers may differ 
significantly due to the orders of magnitude differences in 
their coverage footprints. Therefore, the effective resources Ok 
available in small cells for each user might be significantly 
higher than the macrocells. In such a case, it might be 
beneficial for a user to connect to a small cell even though 
it may not provide the best SIR over the network. We will 
capture this characteristic of HetNets in our definition of rate 
distribution below. Due to the interpretation of a minimum 
rate required by each application, e.g., video, we will study 
rate distribution in terms of "rate coverage", which is defined 
below. It is just the CCDF of rate when the target rates are 
the same for all the tiers. 

Definition 3 (Rate Coverage). A typical user in an open access 
network is said to be in rate coverage if its effective downlink 
rate from at least one of the BSs in the network is higher 
than the corresponding target. We denote the target rate for a 
k th tier BS as Tj.. Rate coverage can now be mathematically 
expressed as: 



Rc — 



( I) max0 fc log 2 (l + SIR(x fc )) > T k ) . (17) 
\keic / 

For closed access, the expression remains the same except that 
the union is now over the set B C /C. 

It is easy to show that the rate coverage for open access 
networks is always higher than the closed access networks. 
This follows from the same arguments that were used in case 
of coverage probability earlier in this section. We now state 
the main ordering result for rate coverage in the following 
Theorem. 

Theorem 2 (Ordering result for rate coverage). For two K- 
tier HetNets with the same resource allocation per user Ok 
for each corresponding tier, the one with system parameters 
{Afc} and {^k} has equal or higher rate coverage than the 
one with system parameters {A^} and {^ k } tf^k < ^'kf or 
all k € /C, and A& > A' k for all k G K, in open access and 
all k £ B in closed access. 

Proof: The rate coverage can be expressed as: 
R c = El ( M max O k log 2 (l + SIR(z fe )) > T k ) , (18) 

where SIR(xfc) can be expressed as: 

Pkhkx k \\xk\\~ a 



SIR(xk) = 



E E Pj9j V \\v\\- a ' 

Pk\\x k \\- a 

E E PjZjk\\y\\- a ' 



(19) 
(20) 



where as in the proof of Theorem [T| we define Zjk = 
9jy/hkx k as me ratio of the two Gamma random variables. 
Now note that the indicator function inside the expectation 
of ( fT8j ) is an element wise decreasing function of Zjk, from 



7 



which the result follows on the same lines as the proof of 
Theorem [T] ■ 

Remark 4 (Same ordering for coverage and rate per user). 
From Theorems^and^ we note that the ordering conditions 
for rate per user are the same as that of coverage probabil- 
ity. Therefore, all the conclusions, including the ordering of 
SDMA, SU-BF and SISO transmission techniques, derived for 
coverage probability carry over to the rate rate per user case 
as well. 

Although it is clear from the above discussion that both SU- 
BF and SISO outperform SDMA both in terms of coverage 
probability and average rate per user, it is important to note 
that we have not yet accounted for the fact that SDMA serves 
more users than both SISO and SU-BF, and may result in a 
higher sum-data rate. To address this, we compare the three 
transmission techniques in terms of ASE in the next section. 



IV. Coverage Probability and ASE Performance 

This is the second main technical section of the paper where 
we derive an upper bound on the coverage probability of 
a typical user in a i^-tier HetNet, where the transmission 
techniques adopted by each tier are characterized in terms 
of the shape parameters A k and ^ k of the Gamma rvs. 
Recall that the coverage probability is formally defined in 
Definition [T] For this analysis, we assume that BS locations 
of each tier are drawn from an independent PPP Q k with 
density X k . Although this model is not as general as the 
one considered in the previous section, it is likely accurate 
in modeling the opportunistic deployment of small cells and 
has been validated for planned tiers, such as macrocells, by 
empirical observations |32| and theoretical arguments under 
sufficient channel randomness (33). The downlink analysis is 
performed at a typical used located at the origin, which is 
made possible by the Slivnyak's theorem (34). Recall that the 
original setup requires that all the BSs of a particular tier are 
either in closed or open access. It is important to note that 
under PPP assumption, the proposed model is applicable even 
if a given fraction of BSs of a particular tier independently 
operate in open or closed access. In such a case, we can divide 
the original tier into two tiers with appropriate densities, which 
is enabled by the fact that independently thinning a PPP leads 
to two independent PPPs. For example, if fraction / of k th tier 
BSs operate in open access and the remaining (1 — /) fraction 
independently in closed access, the k th tier can be divided into 
two tiers modeled as independent PPPs $ fc <i) (open access tier) 
and $^.(2) (closed access tier) with densities A fc (i> = fXk and 
X k (2) = (1 — ,f)X k , respectively. 



A. Upper Bound on Coverage Probability 

Before deriving the upper bound, we first derive an expres- 
sion for the Laplace transform of interference. The result is 
given in Lemma [5] and the proof is given in the Appendix [B] 

Lemma 5. The Laplace transform of interference Ci(s) = 



E [e- sI ], where I = £ fceJC P^yWyt" is 

d(s) = exp -s« E X 3 P f C ( a i *j 
\ ieic 



(21) 



where 



^)=~e(!:>(^ 



a A — ' \ m 

rn — 1 



2 2 

m H , rn 

a a 



(22) 



and B(x,y) = f Q t x 1 (1 — t) y 1 dt is Euler's Beta function. 

Using this result, we now derive an upper bound on the 
coverage probability and the result is given in Theorem [3] 

Theorem 3. The coverage probability of a typical user in a 
K-tier open access HetNet is upper bounded by 



Pc < E XkAk 



keK 



where s Xk = (3 k \\x k \\ a P k 1 and 



E 4\ 



i=0 



(23) 



)da*. (24) 



The upper bound for closed access networks is the same except 
that the summation in (|23[) is over the set B instead of K. 



Proof: We prove the result for open access networks and 
will highlight exactly where the proof will differ for closed 
access networks. Starting with the definition of the coverage 
probability: 



E 



(a) 

< E 



1 ( |J (J SIR(x fe ) > p k 



J2 1 ( sir m > fa) 



E l { p khk xk \\xk\\- a > PkI Xk ) 



.x k e<b k 



(25) 
(26) 
(27) 



where (a) follows from the union bound and I Xk is the 
interference received by the typical user when it is connected 
to the k th tier BS located at x k , i.e., 

7 ^=E E p Mv\\v\\- a - (28) 

jeK ye<ij\x k 

Note that for closed access, the summation in ( |26| and ( pTj ) 
will be over B instead of JC. This is the only difference in 
the proofs of open and closed access cases. Continuing with 
the proof of open access case, since the channel power of 
the direct link is independent of everything else, we can take 
the expectation w.r.t. h kXk inside ( |27] i to write the coverage 
probability as: 



Pc < E E 



keK 



E V(h kXk >p k I Xk \\x k \\ a P k 1 ) 



.x k e® k 



(29) 
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Now we first evaluate the probability P(hkx > z ) as follows: 

, fc>2) Ba^£)s e -. 4 g 1 i (30 , 

where (a) follows from h kx ~ r(Afe, 1), and Y{A k ,z) in 
the numerator is the upper incomplete Gamma function given 



by T(A k ,z) = 



,A k -l„- 



e u du, (b) follows by specializing 



the expression of incomplete Gamma function for the case 
when Afc is an integer. Now denote /3fe||xfc|| Q P,7 1 by s Xk and 
substitute ([30]) in d29]l to get: 



A k -1 

Pc < J2 E E e-^'k J2 

k£K x k £<!> k i=0 



j s x k Ix k ) 



('-0 



A fc -1 



{ s x k Ix k Y 



(31) 



dx k (32) 



X k i 

A fc -1 



= E A feE^i / Ei« k e- s ^(sx k IxJdx k A33) 



keK 



x k m 2 



where (a) follows from Campbell-Mecke Theorem [34|. Now 
note that if A k were 1, the expectation term is just Ci x {s Xk ), 
i.e., the Laplace transform of interference evaluated at s Xk . For 
Aj; > 1, we evaluate the expectation in terms of the derivative 
of the Laplace transform as follows: 



[e-*** 1 ** (sxjxj] ( = ] ai k £{Pf Imi (t)}(s Xk ) (34) 



(>>) 



(SxJ 



<*0*J 



£i* k (s Xk ), 



(35) 



where (a) follows from the definition of the Laplace 
transform and (b) follows from the identity t n f(t) < — > 
( _1 )™5^ /: {/(*)}( s )- Substituting this in §55}, we can ex- 
press the upper bound on coverage probability in terms of 
Laplace transform of interference as follows: 

A fc -1 



Pc<E A * E \ J (-*->)' 

1=0 x k m 2 



S(sx k Y 



£i* k (s Xk )Ax k , 



(36) 



which completes the proof. ■ 
We note that the above upper bound involves a derivative 
of Laplace transform, which makes its numerical evaluation 
difficult. However, it is possible to reduce the upper bound 
to a simple closed form for full SDMA and easy to evaluate 
numerical expressions in the other cases. The simplified result 
is given in the following Corollary. 



Corollary 5. For A& = 1, A k can be reduced to 



A k 

and for A k > 1 to 



2_ 2_ 



(37) 



E ,-i E 



i=0 



(-Sx k fe 



i„-Csi 




Proof: For A& = 1, 

A k = f Cr (/yzfcpp- 1 )^ 

Jx k £R 2 
(») / I own. „a 



dx k . (38) 



(39) 



exp 



(-^Ikfcll 



P k "C \dx k , (40) 



x k & 



where (a) follows from the Laplace transform expression 

derived in Lemma [5] Recall that C = J2jeK ^]PfC(u, ^). 
The closed form expression now follows directly by converting 
the integral from Cartesian to polar coordinates. 

For Afc > 1, using the Laplace transform expression and 
calculating its derivative using Faa di Bruno's formula for the 
composite function (f°g)(s Xk ), with f(s Xk ) = exp (s Xk ), and 
g(s Xk ) = —Cs£ k , the result follows. ■ 

We note that the upper bound is in closed form if Afc = 1 
for all tiers. The result is given in the following corollary. Even 
for Afc > 1, the upper bound can be numerically computed 
fairly easily, especially for small values of Afc. 

Corollary 6. The coverage probability in a K-tier open access 
HetNet with each k th tier BS performing full SDMA to serve 
M k users, i.e., A k = 1 V k G K,, is given by 



Pc < 7T- 



Z)fcg/C ^P k P k ° 



(41) 



For the closed access case, the summation in the numerator 
is over B instead of K,. 

We now comment on the tightness of the coverage proba- 
bility upper bound in various regimes and for various trans- 
mission techniques. 

B. Tightness of the Upper Bound 

For conciseness, we will focus on the open access networks, 
with the understanding that all the arguments remain the same 
for closed access case. Since the bound is derived by using the 
union bound in p6} , the tightness depends upon the number 
of candidate BSs that provide SIR greater than the target SIR. 
Denote this random variable by X({A k },{^ k }), which can 
be expressed as 



X({A k },{^ k }) = Y, E l(SIR(a*)>& 



(42) 



ki£Kx k e<S> k 



The bound holds with equality if there is strictly one candidate 
serving BS for a typical user, i.e., V(X({A k }, {<J>fc}) > 1) = 
0. This is the case in single-antenna HetNets for (3 k > 1, Vfc, 
as shown in [91. For any other general system configura- 
tion, the tightness of the bound depends upon whether the 
probability ¥(X({A k }, {^ k }) > 1) is close to zero or not. 
In general, it is hard to evaluate simple expressions for this 
probability. However, it is possible to make a few simple 
observations about the expected tightness of the bound. For 
instance, the bound gets tight with the increasing values 
of target SIRs because X({A k }, {^fc}) is an element-wise 
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decreasing function of /3k- For further insights, we derive 
the following ordering result for X({Ak}, {^>k})- The proof 
follows using Lemma[3]on the same lines as that of Theorem[T] 

Theorem 4 (Ordering result for X). If A k > A' k and 

*fc < *' fe V k, then X({A k }, {%}) (first order) stochastically 
dominates X ({A' k } , {%}), i.e. ¥(X({A k },{^k}) > n) > 
¥(X({A' k },{%})>n),yn. 

Proof: Using the alternate expression of SIR given by |6|, 
express X({A k },{^> k }) in terms of the channel power gains 



E E i 



Pk\\x k \\- a 



>Pk), (43) 



where Zjk = gj y /hk Xk is the ratio of the two Gamma 
random variables. For another system with Z'^ k = g'j y /h' kXk , 

Zjk <st Z jk if ^fe ^ ^k and *j - which follows from 
Lemma |3] The result now follows on the same lines as the 
proof of Theorem [TJ using Lemma [4] along with the fact that 
X({ Ak}, {^k}) is an element-wise non-increasing function 
of Z jk . ■ 

Remark 5 (Tight Bound in case of SDMA). One of the useful 
consequences of Theorem^is the prediction of the tightness of 
the upper bound for SDMA. One interpretation of the above 
result is that the bound gets tighter when all the BSs serve 
more users, i.e., Ak decreases and ^!k increases for all the 
tiers. A limiting case is that of full SDMA, where the number of 
users served by each BS is equal to the number of its transmit 
antennas. Beyond this point, the bound gets tighter with the 
addition of more transmit antennas keeping Ak = 1. We revisit 
these observations in the numerical results section and show 
that the bound is in fact surprisingly tight even for two transmit 
antennas down to very low target SIRs. 



In the rest of this section, we will mainly focus on the full 
SDMA case. Recall that in this case A k — 1 and ^ = Mj and 
the coverage probability upper bound is given by Corollary [6] 
As argued in Remark [5] and validated in the numerical results 
section, the closed form upper bound is tight and can be 
used as an approximation for the coverage probability. For 
simplicity we will use equality instead of an approximation. 

Remark 6 (Similarity with P c in SISO case). The coverage 
probability expression derived for full SDMA case in Corollary 
[6] has a striking similarity with the coverage probability in 
the single-antenna case derived in /[9^. The only difference is 
that the constant C(a,Mj) in that case is simply C(a) = 

2ir 2 csc(^) 
a 

To facilitate direct comparison of the full SDMA and the 
SISO cases, we need to understand the relationship between 
C(a) and C(a, M). Let us take a closer look at the expression 
of C(a,M) given by ((22]). First note that C(a,M) is an 
increasing function of M. Now let us evaluate C(a, 1): 



C(a,l) = — b( 2 ,1- 2 
a \ a a 



2tt /2\ / 2\ 2tt 2 csc(^) 

— r _ r i - - = ^a, (45) 

a \ a J \ a J a 



where the last step follows by Euler's reflection formula. 
Hence C(a, 1) is the same as C(a) derived for the single- 
antenna case in 19). From the monotonicity of C(a,M) it 
follows that C(a,M) > C(a) VM > 1. 

Remark 7 (Full SDMA vs. SISO coverage). Keeping all the 
system parameters the same, the full SDMA coverage is always 
lower than that of the SISO case. This is consistent with the 
coverage probability ordering results derived in the previous 
section. 

Remark 8 (Scale invariance in open access HetNets). The 
full SDMA coverage probability is invariant to the density 
of the BSs, number of tiers and the transmit powers when 
the target SIRs and the number of transmit antennas are the 
same for all the tiers in open access HetNets. The coverage 

2 

probability in this case is given by P c = c{a M) @~™ • 
result is again similar to the single-antenna result where the 



coverage probability reduces to P c 



The scale 



invariance result does not hold for closed access HetNets. 



C. Area Spectral Efficiency 

Although the comparison of various system configurations 
and transmission techniques is quite conclusive in terms of 
coverage probability and the rate per user, it does not directly 
account for the fact that some techniques, such as SDMA, 
serve higher number of users than the others, such as SU-BF, 
and may result in higher sum data rate. To account for this fact, 
we consider ASE, which gives the number of bits transmitted 
per unit area per unit time per unit bandwidth. For a multi-tier 
setup, it can be formally defined as 



V= ^*fc A fc l0 g2(l+A)P< 



(k) 



(46) 



where P c ^ fe ^ is the per tier coverage probability, i.e., coverage 
probability conditional on the serving BS being in the k th tier. 
Since the derivations of per tier coverage probabilities are out 
of the scope of this paper, for analytical comparisons we limit 
our discussion to the cases where — P c for all tiers. 

This is guaranteed for the SISO case when the target SIRs 
are the same for all tiers and for SDMA when additionally 
the number of antennas per BS are also the same for all tiers. 
Recall that the coverage probabilities under these assumptions 
are scale invariant, as discussed in Remark [8] The ASE under 
these assumptions can be expressed as: 



7 ? = P C log 2 (l + /3)^* fc Afc. 



(47) 



keic 



We first compare the ASE of the full SDMA and the SISO 
cases below. The ASE for full SDMA case is: 

VM = M * /T° log 2 (l +p) X k, (48) 

' k£K 



(44) and for the single-antenna case is: 



Vs 



C(a) 



p-- bg 2 (l +/3)E 

k£K 



(49) 
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The ratio of the ASEs can be expressed as: 

77m = MC(a) 
r] S C(a,M)' 

Using the fact that 

C(a,M) . , 

lim v ' ' = ttT 1 - 2 a), 

M-coo Mi 

the ratio of the ASEs can be approximated as: 

r, M „ M^iCja) 
Vs ~ 7rr(l - 2/a) 

= v(l + -]M 1 -i, 



(50) 

(51) 

(52) 
(53) 



which shows that the ratio grows with the number of antennas 
when a > 2. In the next section we will validate this 
observation and show that the ASE in case of full SDMA case 
is always higher than the SISO case. We will also show that 
the approximation is surprisingly tight even for small values 
of M. 

Another relevant comparison is that of full SDMA and SISO 
transmissions when both the systems are serving the same 
density of users. To facilitate this comparison, the densities 
of BSs for SDMA case will be lower than the SISO case 
by a factor of M. This comparison will provide insights into 
whether it beneficial in terms of ASE to deploy A BSs per unit 
area with M antennas or MA single-antenna BSs per unit area. 
In that case, the ratio of the ASEs can be approximated as: 



«-SDMA: 
e - SDMA: 
SDMA: 
$ - SDMA: 



M=2 Simulation 
M=2 Theory 
M=4 Simulation 
M=4 Theory 



is 



r i 



2 |M- 



(54) 



which shows that the ratio decreases sublinearly with the 
number of antennas when a > 2. We will validate this 
observation in the next section and show that the ASE in case 
of SISO is higher than the full SDMA case. 

So far, we have focused only on the comparison between 
full SDMA and SISO cases, mainly because the coverage 
probability expressions for these cases are known in closed 
form. Since the coverage probability upper bound for SU-BF 
cannot be reduced to closed form and moreover the tightness 
of the bound is questionable, we cannot perform similar 
comparisons with SU-BF unless we derive a simple coverage 
probability expression, which is out of the scope of this paper. 
Having said that, it is possible to compare the three cases in 
the very low target SIR regime, where the coverage probability 
approaches 1. Therefore, the ASE of SISO and SU-BF cases 
should be similar because they serve the same number of users 
and of SDMA case should be higher than the both when the 
densities of BSs are the same for all the cases. We will revisit 
this case along with the moderate and high target SIR regimes 
in the next section. 

V. Numerical Results 

Since there is a slight difference in the simulation of the 
proposed multi-antenna model and the ones proposed in the 
literature for the single-antenna HetNets, e.g., 15), we will 
briefly summarize the simulation procedure before explaining 
the results. Choose a sufficiently large window and simulate 
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Fig. 3. Coverage probability of a two tier HetNet when both the tiers perform 
full SDMA (K = 2,P = [1..01], Mi = M 2 = M,X 2 = 2Ai,ft = 
/8 a , a = 3.8). 
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Fig. 4. The ratio of the full SDMA and the SISO ASEs given by J50J. The 
dotted plots correspond to the approximation of the ratio given by {52} . 



the locations of different classes of BSs as realizations of in- 
dependent PPPs of given densities. Associate two independent 
marks h x and g x with each BS. Assuming the typical user lies 
at the origin, calculate the desired signal strength from each 
BS using the sequence of marks {hk x } and the interference 
power using the sequence {gkx}- Calculate the received SIR 
from each BS. The user is now said to be in coverage if 
the received SIR from at least one of the BSs belonging to 
the permissible tiers is more than the corresponding target. 
Repeating this procedure sufficient number of times, we have 
an estimate of the coverage probability. For conciseness, we 
focus only on the open access case in this section. 

Using this procedure, we now evaluate the coverage proba- 
bility of a two tier HetNet in full SDMA case and compare the 
results with the upper bound derived in Corollary [6] in Fig. [3] 
As stated in Remark [5] the bound is tight down to very low 
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Fig. 5. Comparison of the ASE in a two tier HetNet (K = 2, P = 
[1,.01],A2 = 2Ai,/3i = (52, a = 3.8). The number of antennas in case 
of multi-antenna tiers is M = 4. Full SDMA corresponds to \I/ = M. 



Fig. 6. Comparison of the coverage probability in a two tier HetNet for var- 
ious combinations of multi-antenna techniques (K = 2, P = [1, .01], A2 = 
2Ai,/3i = p2,a = 3.8). The number of antennas in case of multi-antenna 
tiers is M = 4. 



target SIRs. Even for M = 2, the bound is tight down to about 
—4 dB. A slight gap at moderate to high target SIRs is due 
to the edge effects in simulation, also observed earlier in [9|. 
This validates our assumption of considering the upper bound 
as an approximation of the coverage probability in case of full 
SDMA in the previous section. 

We next plot the ratio of the ASEs of the full SDMA and 
SISO cases given by |50} in Fig. [4] As stated in the previous 
section, the ASE of full SDMA case is always higher and the 
gap increases further at higher path loss exponents or when 
the number of transmit antennas increases. We also note that 
the approximation given by ( f52] > is extremely tight even when 
the number of antennas is small. 

We now additionally consider SU-BF and compare ASE for 
the three transmission techniques in a 2-tier HetNet in Fig. [5] 
Both the tiers are assumed to follow the same transmission 
technique and the ASE result for SU-BF is computed nu- 
merically by computing the per tier coverage probability. For 
comparison, we consider two cases, one in which the density 
of the BSs in the three setups remain the same, and the other in 
which the densities are adjusted such that the density of users 
served in the three cases is the same. In the first case, SU- 
BF, which always outperforms SISO, even outperforms full 
SDMA in the high target SIR regime despite serving smaller 
number of users. The trends in the low target SIR regime are 
consistent with arguments presented in the previous section. In 
the second case, where the density of the users is the same in 
all the cases, the ordering of the three transmission techniques 
in terms of ASE is the same as that of coverage and rate per 
user. 

Before concluding this section, we study the effect of adding 
a second tier to the network in terms of coverage probability in 
Fig.|6j where both the first and the second tier can be one of the 
three possible types: i) SISO, ii) full SDMA, iii) SU-BF. This 
result shows that the case where both the tiers perform SU-BF 
results in the highest coverage, whereas the case where both 



the tiers perform full SDMA leads to the lowest coverage. This 
is because SU-BF case has an additional beamforming gain; in 
addition to the proximity gain enjoyed by the SISO case. This 
result is consistent with the coverage ordering results derived 
in this paper. 

VI. Conclusions 

In this paper, we have proposed a new tractable downlink 
model for multi-antenna HetNets. For any given BS distri- 
bution, we derived ordering results for coverage probability 
and per user rate to compare different transmission techniques, 
such as SDMA, SU-BF and baseline single antenna transmis- 
sion. In addition to significantly generalizing the state of the art 
PPP based random spatial models for cellular networks, this 
approach circumvents the need for deriving explicit expres- 
sions for coverage and rate, which may not reduce to simple 
closed forms in all the cases. Our analysis demonstrates that 
for a given total number of transmit antennas, it is preferable 
to spread them across many single-antenna BSs vs. fewer 
multi-antenna BSs, both in terms of coverage and rate per 
user. We also show that SU-BF provides higher coverage 
and rate per user than both SISO and SDMA due to an 
additional beamforming gain. To account for the fact that 
certain transmission techniques, such as SDMA, serve more 
users and may provide higher sum-rate, we derive an upper 
bound on the coverage probability assuming an independent 
PPP model for BS locations and use it to compare different 
transmission techniques in terms of ASE. 

Future work could consider HetNets with multi-antenna re- 
ceivers and investigate potential gains by performing interfer- 
ence cancelation and/or receive combining. Further extensions 
to this work could include the effect of imperfect CSI at 
the transmitters and opportunistic scheduling on the coverage 
probability and spatial reuse. Another important extension of 
the modeling tools developed in this paper is to the uplink of 
cellular networks (35). 
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Appendix A 
Signaling Preliminaries 

The received signal y k from k th tier BS at a typical user 
located at the origin is given by 



Vk 



2 vl 



E E 

keKy£<S> j \x k 



(55) 



where P k is the per user transmit power of k tier BS, and 
is the normalized transmit signal vector. The 



Zfe G 



■>Af k Xl 



' (V fc V fe " )- x e C MfcXWfc , for V = 
yS k xM h jjgjjjg fhe concatenated 



channel vector from k th tier BS to a typical user located at 
origin is denoted by Vkx G C A/fcXl and for the interfering 
link from a j t/l tier BS located at y E M. 2 is denoted by uj y £ 
C Mj Xl . The vectors v, u are assumed to have i.i.d. CAf(0, 1) 
entries, independent across BSs and of the user distances. 

This paper assumes linear precoding, in which the k th 
tier BS multiplies the data symbol s k .i destined for the i th 
user, for 1 < i < by Wk,i so that the transmitted 

signal is a linear function, i.e. z k = ^k.iS k ,i- When zero- 

i=l 

forcing beamforming with perfect CSI is employed to serve 
\E , / C , tyj users in tier k, j respectively, the columns of the 
precoding matrix = [wk,i]i<i<y k G C MfcX * fc equal the 
normalized columns of Vfe (^r.\r. , i .- ,r 1 w... >j<,, 
[v 1; ...,v fc ... v*J* G 
matrix of channel directions, where the direction of each vector 
channel is represented as v^. = The desired channel 

power is given by h kx = |v^w fcjfc | 2 = \v* kx Mv ktk \ 2 ■ \\v kx \\ 2 
which equals the product of two independent rv's which are 
distributed as Beta(Af fc - * fc + l,^ k - 1) and T(M k ,l), 
respectively. Therefore,the channel power is h kx ~ T(Afc, 1) 
with A k = Mk — \?k + L For the distribution of the interfering 
marks, we assumed that the precoding matrices have unit- 
norm columns and that Wj is calculated independently of 
\ij y . Therefore, iij y and are independent isotropic unit- 
norm random vectors, and |u* y Wj | is a linear combination 
of vffj complex normal random variables, i.e. exponentially 
distributed. As a result, we have that gj y ~ r(^J/j, 1), since it 
is the sum of fyj i.i.d. exponential random variables. 

The case A k = 1 and ^ = Mj is referred to as full SDMA. 
The case that each BS only serves one user, i.e. ^> k = 1, using 
the beamforming vector w kx = v kx corresponds to SU-BF or 
MISO eigen-beamforming. In that case, the channel power is 
given by h kx ~ r(Afc, 1) with A k = M k and the interference 
marks as gj y ~ r(\&,-,l) with fyj = 1, Vj G IC, since the 
beamforming vectors Wj y used by the j th tier interfering BS 
are calculated based on vj, i.e. independently of Uj y . 

Appendix B 
Proof of Lemma[5] 

The Laplace transform of interference £/(s) = E/ xfc [e -57 ] 
can be derived as follows: 



E 

(a) 



n 



-*Pj9jv\\v\\~ 



(56) 
(57) 



n Kibir a 

= } [] cxp f-A j / (1 - ( S Pilb||- a )) dy) (59) 
j&c V • /r2 ' 




n ex p 

3&C 



(60) 

(i+ s p,n y ii-)^ dy ; (oi) 

(IT^NF^ dy J (62) 

(63) 
(64) 



(Z) JJ e -*rAi(«- fc ^)= E m ta (*>)/a°° rdr 
J'6/C 



(9) 



je/c 



(65) 



where (a) follows from the independence of the tiers, (6) 
follows from the fact that channel powers are independent 
of the BS locations, (c) follows from PGFL of PPP [34], (d) 
follows from the Laplace transform of the gj y ~ r(4 , j,l), 
(e) follows from Binomial theorem, and (/) follows from 
converting to Cartesian to polar coordinates, and (g) follows 
by substituting (1 + r - ") -1 — > t to convert the integral into 
Euler's Beta function B(x, y) = t x - x {l - t) y - x dt. 
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